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With the help of the replica exchange Monte 
Carlo method and the improved Monte Carlo 
renormalization-group scheme, we investigate 
over a wide area in the phase diagram of the 
Gaussian random field Ising model on the sim- 
ple cubic lattice. We found that the phase tran- 
sition at a weak random field belongs to the same 
universality class as the zero-temperature phase 
transition. We also present a possible scenario 
for the replica symmetry-breaking transition. 
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The random field Ising model in three 
dimensionB has been intensively studied for 
over 30 years by experimental, theoretical, 
and numerical methods. But there are still 
some unanswered questions. One of the ques- 
tions is as follows: Does the phase transi- 
tion for different strengths of random field 
belong to the same universality class? Crit- 
ical behavior at the strong field region!!! 
and zero-temperature lineS has been stud- 
ied numerically by several authors, and from 
the fact that the observed value of the crit- 
ical exponent 9 is positive]! the existence 
of the zero-temperature random fixed point 
(ZRFP) seems unquestionable. Together 
with the fact that random field is relevant at 
the pure critical point ,0 the renormalization- 
group (RG) flow in the temperature-field 
phase diagram can be depicted as in Fig. |l|. 
The behavior of the RG flow in the inter- 
mediate region is unknown and it is unclear 
whether the ZRFP governs the whole phase 
transition or not. 



Another question is as follows: Is there a 
glasslike phase near the critical line in the 
phase diagram? Theoretical studies based 
on the replica formalism! have predicted the 
presence of a glasslike phase which is char- 
acterized by replica symmetry breaking near 
the transition temperature.! But owing to the 
unphysical n — > limit in the theory, inter- 
pretation of the result is not trivial. 

For the investigation of these problems, 
direct observation of the RG flow may be 
the most powerful method. In the present 
work, we numerically observed the RG flow 
of the Gaussian random field Ising model over 
a wide region in the phase diagram using 
an improved Monte Carlo renormalization- 
group (MCRG) scheme,! which is a very sim- 
ple and sophisticated method. We used the 
following model in the Monte Carlo simula- 
tions: 
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where the first summation runs over all the 
nearest-neighbor pairs of the sites on an L x 
LxL simple cubic lattice (up to L = 16) with 
a periodic boundary condition, and hi is a 
random Gaussian variable whose average is 
and variance is 1. We measured the following 
quantities: 
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where M denotes magnetization, the sym- 
bol (• • ■) and [■ • ■] denote thermal and en- 
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semble averages, respectively. Tl and Sl 
are the amplitudes of thermal and sample- 
to-sample fluctuations, respectively, of the 
rescaled block spin Mj J[(M 2 )] defined on an 
L x L x L block.ll Thus Tl and Sl reflects 
the renormalized temperature and strength 
of the random field, respectively. If the ZRFP 
governs the whole critical line, the flow of 
T L and Sl will become the one depicted in 
Fig. 0(a). On the other hand, if there ex- 
ists another fixed point, the flow will become 
the one depicted in Fig. @(b). Note that, in 
the paramagnetic phase, thermal and sample- 
to-sample fluctuations are on the same order 
0(L 3/2 ) and there are a set of "high temper- 
ature fixed points" on a line T L + Sl = 1, 
where both fluctuations are Gaussian. The 
position of the fixed point depends on the ra- 
tio between the two amplitudes. 

In the Monte Carlo simulation, we used 
a single-spin-flip Metropolis update and the 
replica exchange method.! In a replica ex- 
change procedure, an exchange for all adja- 
cent temperature pairs was tried. The ex- 
change procedure was performed for every 
two Metropolis sweeps. We used 20 replicas 
for all field realizations. Temperatures for all 
replicas were chosen so that the replica ex- 
change occurs with modest frequency and the 
highest temperature is well above the tran- 
sition temperature where the single-spin-flip 
provides good ergodicity. Figure |3] shows the 
region in the phase diagram where the simu- 
lations were carried out in the present work, 
together with other recent numerical stud- 
ies. The case h = 1.3 was the hardest to 
relax, in which we used 88 000 Monte Carlo 
steps (discarding the initial 8000 steps) for 
each field realization. We checked that two 
kinds of initial states, all spins up and down, 
give the same thermal averages of M within 
statistical errors, thus the system is well equi- 
librated. This result is remarkably fast com- 
pared to the simulation in which only the 
Metropolis update was used and 2 x 10 6 steps 



were needed to fulfill the same criterion.! 2 ] We 
used 120 different random field realizations 
for each L and h. This number of samples 
is rather small compared to other recent nu- 
merical works. But our aim in the present 
work is to draw qualitative conclusions, not 
quantitative estimates of critical exponents, 
and it was enough for this purpose. 

Figure [| shows the flow of T L and Sl on 
a logarithmic scale. All lines are drawn from 
(T L , S L ) to (T 2 l, S 2 l) with L = 4 (thin lines) 
and L = 8 (bold lines). The RG flow of both 
sizes are consistent, which indicates that the 
finite-size effect is small enough and the final 
destination of the RG flow can be safely pre- 
dicted. Thus one can see that there are no 
fixed points in the intermediate region, and 
the renormalization flow which starts at the 
weak-field region is eventually attracted to 
the low-temperature region T L < 0.01. Fig- 
ure [5] shows the RG flow near the ZRFP. 
The precise position of the ZRFP (0, S 1 *) is 
hard to estimate, owing to its proximity to 
the low-temperature fixed point (0,0). @ We 
only performed a rough estimate of the criti- 
cal value S 1 * as < 5** < 0.01, which is consis- 
tent with the result of the zero-temperature 
simulation,!^ = 0.003(2). Thus Figs. | and 
[5] indicate that the ZRFP governs the whole 
phase-transition line. 

Now let us consider the possibility of a 
first-order transition. If the transition is of 
first order, an ensemble distribution of (M 2 ) 
will exhibit a double-peak structure and a 
plot of Sl will develop a sharp peak at the 
transition temperature. But no such behav- 
ior was observed in the present work. Note 
that when we observe a single sample at 
zero temperature and varying h, the sys- 
tem jumps between different ground states 
and there are many discontinuous transition 
points on the h axis. Similar behavior will 
be observed in the low-temperature region 
and somes sample may exhibit the double- 
peak behavior.S However, if ensemble distri- 
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bution of the ground-state magnetization at 
zero temperature does not exhibit double- 
peak behavior, an ensemble average of the 
ground state magnetization will not change 
discontinuously and the transition is of sec- 
ond order, as is the result of RefM 

Figure |^ shows plots of Sl against tem- 
perature at h — 1.3. The bold line denotes 
the expected behavior in the L — > oo limit, 
assuming that the basin of the fixed point 
Sl = 1,Tl = has a finite measure in 
Fig. 0(a). Soo reaches its maximum value 
1 slightly above the critical temperature and 
remains 1, then discontinuously drops to a 
critical value S* at the critical temperature. 
Note that Sl is a similar quantity to the 
replica interaction term 2 0| in the 4 model 
based on the replica formalism,! in which 
divergence of the coefficient of a ^ the 

paramagnetic phase is reported.! In the S^ = 
1 region, thermal fluctuation of the magne- 
tization becomes more and more negligible 
compared to its absolute value as L becomes 
large. This means that the breadth of each 
valley in the free-energy landscape becomes 
negligible compared to the distances between 
each valley. In this limit, an infinitesimal 
change of h or T will drive the system out 
of a valley and into another valley.0 This 
could be a possible scenario for the replica 
symmetry-breaking transition!!!'! and can be 
tested numerically by checking the following 
inequality: 



lim lim \{M(h))(M(h + A))] 

* i(M(h)n 
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with a very careful finite-size scaling analysis. 

In conclusion, we performed a replica ex- 
change Monte Carlo simulation of the Gaus- 
sian random field Ising model in three dimen- 
sions and, using an improved MCRG scheme, 
we found that the phase transition at the 
weak random field belongs to the same uni- 
versality as the zero-temperature phase tran- 
sition, and we presented a possible scenario 



for the replica symmetry-breaking transition. 
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FIG. 1. Schematic renormalization flow in the 
temperature-field phase diagram. 
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FIG. 2. Expected renormalization flow of Tl 
and Sl for the case of (a) the zero-temperature 
random fixed point only and (b) another fixed 
point in the intermediate region. L, P, and R x 
denote low-temperature, pure-critical, and ran- 
dom fixed points, respectively. 
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FIG. 3. Regions in the T-h phase diagram 
where simulations were carried out in the present 
work, together with other recent numerical 
works. 
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FIG. 4. Observed renormalization flow of Tl 
and Sl- The bold dashed line is an approximate 
phase boundary. Note the logarithmic scale of 
the plot. Errors are smaller than the size of sym- 
bols, unless errorbars are explicitly shown. 
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FIG. 5. Renormalization flow of Tl and Sl 
near the ZRFP. The bold dashed line is an esti- 
mated position of ZRFT in Ref.[4]. 
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FIG. 6. Plots of Sl against temperature T for 
h = 1.3. The bold line is the expected behavior 
for the L = oo case. 
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